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Abstract. Recent experiments have shown that it is possible to study a fundamental astrophysical process such as 
dynamo action in controlled laboratory conditions using simple MHD flows. In this paper we explore the possibility 
that Taylor-Couette flow, already proposed as a model of the magneto-rotational instability of accretion discs, can 
sustain generation of magnetic field. Firstly, by solving the kinematic dynamo problem, we identify the region of 
parameter space where the magnetic field's growth rate is higher. Secondly, by solving simultaneously the coupled 
nonlinear equations which govern velocity field and magnetic field, we find a self-consistent nonlinearly saturated 
dynamo. 
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1. Motivation 

Even though astrophysical objects possessing observable 
magnetic fields are extremely diverse, it is widely accepted 
that the physical mechanisms supporting the fields are 
fairly universal and rely on features common to virtually 
all astrophysical objects (e.g. differential rotation, con- 
vective or turbulent motions, etc.). MHD dynamo theory 
quantifies this idea and states that astrophysical magnetic 
fields are created by inductive currents driven by motions 
of electrically conducting plasmas (Moffatt 1978). Until 
recently, dynamo action was the subject of theoretical or 
numerical investigations only. The recent demonstration 
of dynamo action in controlled laboratory experiments 
(Gailitis et al. 2001; Stieglitz & Muller 2001) has stim- 
ulated interest in the study of dynamo action in confined 
geometries of potential laboratory interest, such as spheres 
or cylinders. The configuration of an MHD fluid confined 
between concentric cylinders (Taylor-Couette flow) is par- 
ticularly relevant. Dynamo action is associated with spiral 
and sheared flows, which suggests that Taylor-vortex flow 
is good candidate for dynamo experiments in a simple ge- 
ometry. 

At the same time, Taylor-Couette flow is already a 
useful benchmark for studying instabilities relevant to as- 
trophysical processes (Balbus & Hawley 1991; Ji et al. 
2001; Riidiger & Zhang 2001), as the flow pattern can 
model Keplerian motion. In particular we note the works 



of Dobler et al. (2002) who studied the dynamo mecha- 
nism in the presence of an imposed axial flow (result- 
ing in a screw dynamo) and of Laure et al. (2000) who 
concentrated on the basic Taylor-Couette configuration. 
The latter performed a kinematic dynamo calculation in 



this geometry and demonstrated that an imposed Taylor- 
vortex flow is capable of dynamo action. The combina- 
tion of shear with simple roll flows has also been mod- 



elled by Dudley & James (1989) in spherical geometry. 



There it was demonstrated that magnetic field genera- 
tion is sensitive to the nature o f the driving flow. Our 
wor k differs from the st udies of Dudley fc James (1989) 
and Laure et al. (2000) in two important respects: first 
(kinematic dynamo) we use the velocity fields that are ac- 
tual solutions of the Navier-Stokes equations to generate a 
magnetic field, not arbitrary imposed flow fields; secondly 
(fully self-consistent dynamos) we let these velocity fields 
evolve alongside the magnetic field, thus solving the full 
MHD equations. 

The a im of this work is tw o-fold. Firstly we widen the 
study of Laure et al. (2000) in the parameter space by 
finding the growth of the magnetic field as a function of 
the flow patterns. Secondly, and more importantly, we go 
beyond the limits of (linear) kinematic theory and inves- 
tigate the fully self-consistent dynamo mechanism (mag- 
netic fields and flow pattern affect each other and saturate 
nonlinearly) . 

2. Model 

We consider an incompressible fluid contained between 
two coaxial cylinders of inner radius R\ and outer ra- 
dius i?2 which rotate at prescribed angular velocities fii 
and 0,2- The velocity and magnetic fields, V(r,9, z,t) and 
B(r,0, z,t) are determined by the MHD equations which 
we write in dimensionless form as 



dV 

dt 



(V ■ V)V = -Vp + V 2 V + (V A B) A B, 



(1) 



2 
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dB 

~dt 



P 



-V 2 B + VA(VAB), 



V-B = 0, V-V = 0, 



(2) 
(3) 



where p is the pressure. In writing (|l|) and (^) we used 
S = R 2 — Ri as unit of length, <5 2 /V as unit of time and 
(flop) z v/8 as unit of magnetic field, where the kinematic 
viscosity u, the density p and the permeability po are con- 
stant. The governing dimensionless parameters of the sys- 
tem are the Reynolds numbers Rei, Re2, the radius ratio 
77, and rotation ratio p, 

RiQxS R 2 Q 2 6 Ri n 2 

Rei = . Re-2 = , ri = —, p=—,(A) 

v v R 2 Hi 

together with the magnetic Prandtl number P m , 
P,„ = - (5) 



v 
A ' 



where A = \j(poo) is the magnetic diffusivity and a is the 
electrical conductivity. 

We assume no-slip boundary conditions for V, and 
electrically insulating boundary conditions for B. For 
modes of the form B(r,t,z) = B(r) e 1 ( Q ^+ m8 ) these con- 
ditions are 



a = m = 
a — 0, m 7^ 

a ^ 



Bg = B z = 0; 

±TO-B# 

flUB, 



B, 



B z = 0; 
r 



(G) 



where for ± we take + at — at R 2 . The symbol B m (r) 
indicates the modified Bessel functions, I m (ar) at Ri and 
K m (ar) at i?2- The axial wavelength of a pair of Taylor- 
vortices is 2ir/a. 

Our numerical method for timestepping the 
3D nonlinear MHD equations is detailed in 

The formulation is based 
B with the toroidal-poloidal 



Willis fc Barenghi (2002) 



on representing V 
decomposition 



and 



A = ip 6 + 4>oZ + VA (ipr) + V A V A (>r), 



(7) 



where tp(r, t, z), <p(r, t, z) and ipo( r )> 4>o( r ) contain the pe- 
riodic and non-periodic parts of the field respectively. The 
potentials are expanded spectrally over Fourier modes in 
the azimuthal and axial directions and over Chebyshev 
polynomials in the radial direction. 

The governing equations for the magnetic field are the 
r-components of the induction equation and its first curl. 
For the velocity we follow the procedure applied to the 
magnetic field and take the r-components of the momen- 
tum equation and its first curl. As the pressure has not 
been eliminated, we also take the divergence to obtain 
the pressure-Poisson equation. All five governing equa- 
tions are second order in r. Time stepping is based on 
a combination of second order accurate Crank-Nicolson 
and Adams-Bashforth methods. The code has been tested 
against published results with and without a magnetic 
field (Chandrasekhar 1961; Roberts 1964; Marcus 1984; 
Jones 1985; Barenghi 1991). 



tj = 0-5 



Fig. 1. Growth rates as for various r\ and m, where Rei = 
2i?e lc (r7), P m = 2. 

3. Kinematic dynamo 

For prescribed Reynolds number Rei , wavenumber a, ro- 
tation ratio p and radius ratio 77 the Navier-Stokes Eq. (Q) 
is solved in the absence of a magnetic field. We always as- 
sume Rei > Reic, where Rei c is the critical Reynolds num- 
ber at which azimuthal Couette flow is unstable to the for- 
mation of axisymmetric Taylor vortices (Chandrasekhar 
1961). We also assume a = 3.14, which corresponds to 
almost square cells. The nonlinear Taylor-vortex flow V 
thus obtained is used when solving the induction Eq. ^ 
for B, starting from a small magnetic seed field of wave- 
length as- 

Note that V being fixed, (||) is linear and has eigen- 
function solutions B which grow or decay exponentially. 
If the real part of the growth rate <jb is positive then 
the magn etic field grows (k inematic dynamo action). 



Following [Laurc ct al. (2000) 



we assume olb = ^a, 



the 



characteristic length for the magnetic field is the length of 
two pairs of Taylor-vortices. In several calculations with 
olb — a we found the magnetic seed field decayed quickly. 
The axisymmetric flow cannot generate an m — mag- 
netic field. In Fig. ^ we see that in narrower gaps the 
dynamo prefers larger m. The dynamo is local here in the 
sense that the characteristic length scale for the magnetic 
field comply with the scale of the flow. As the Taylor- 
vortex flow is itself unstable to non-axisymmetric pertur- 
bations in narrower gaps (see Jones 1985), we usually con- 
sider the case rj = 0.5 where m = 1 is preferred. The 
relative stability of Taylor-vortex flow in wider gaps also 
allows for a clearer interpretation of the affect of the mag- 
netic field in the nonlinear self-consistent solutions pre- 
sented in Sect. || 

Figure || shows that, not surprisingly, ob falls off 
sharply when P m < 1 , in which case much larger Reynolds 
numbers (relative to Rei c ) are needed for dynamo action. 
Generally it is easier to generate a magnetic field in media 
with larger magnetic Prandtl numbers. This is believed to 



be the case for warm weakly ionised gas, Mclvor (1977) 
Kulsrud fc Anderson (1992) 
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Fig. 2. Growth rates a B for various 77 as a function of Fig. 4. Growth rates os with co-rotation, ri = 0.5, P m = 
P m , where itei = 2Rei c (rj), m = 1. Critical Reynolds 2. At = 0, ^t? 2 , §77 2 , |?7 2 , critical Reynolds numbers are 



numbers at 77 = 0.3,0.5,0.7,0.9 are respectively Re\ c 
73.3,68.2,79.5,131.6. 



-1.0 1 

1.0 




Fig. 3. Growth rates 05 as versus Re\ at 77 = 0.5. 



Figure [| shows growth rates as a function of Re\ with 
R&i — 0, 77 = 0.5 for a few values of P m . Driving the 
flow harder seems to move the flow into a regime which 
is less favourable for magnetic field generation when the 
outer cylinder is fixed. Allowing co-rotation, the stabil- 
ity boundary for onset of Taylor-vortex flow tends to the 
Rayleigh line fj, = rf 2 for large Reynolds numbers. As 
circular-Couette flow O(r) = a + b/r 2 alone is not ca- 
pable of dynamo action, we must consider flows which are 
Rayleigh unstable, [i < rj 2 . Figure || shows growth rates as 
a function of Re\ for different values of /i. Again, we nor- 
malise the intensity of the drive Re\ by the critical value 
Re\ c . It is apparent that magnetic field generation is easier 
with co-rotating cylinders and /j = ^if is an important 
case. At this particular (i the amplitude of the azimuthal 
disturbance Vg = Vg — rfl(r) and V z are similar, seen 
in Fig. ||. At a given Re% / Rei c (fj.) , the components V r , 
V z appear to be very similar in structure and amplitude. 
However, Vg changes significantly with \x. 

Of the rotation ratios considered in Fig. ^, /1 = |?7 2 
would appear to be most favourable for magnetic field gen- 



Re lc = 68.2, 84.3, 112.7, 155.3 respectively. 
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5. Amplitude of the azimuthal disturbance Vg — 
r fl(r) drops as /j, — > rj 2 in the absence of a magnetic 
field. At each Re\/ Re\ c the amplitudes of V z are close for 
all (i. 



Fig. 6. Dependence on P n , 
field generation. \i — jrf, rj 



for the driving required for 
= 0.5. 
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Fig. 7. Magnetic energy versus time as the magnetic field 
saturates, r\ = 0.5, P m = 2. (a) fj, = 0, Rei = 1.5i?ei c , 
Eccf = 4.15 x 10 4 . (6) n = \r, 2 , Rei = 2Re lc , E ( 
3.63 x 10 5 . 



'CCF 



eration. We define Re\ as the critical Reynolds number at 
which erg = (marginal state for the onset of dynamo ac- 
tion) for this [j>. Figure |6| shows how Re\ depends on P m . 
Fitting the last three points for small P m we obtain the 
slope —1.1. Approximately, Re[/Re\ c oc 1/P m , for small 
P m . Our interest in this case is linked to the small mag- 
netic Prandtl numbers of laboratory fluids. Defining the 
magnetic Reynolds number Re m as 



Re m = = Rei P„ 



(8) 



we conclude that the critical magnetic Reynolds number 
Re^ = Rei P m is approximately constant and O(10 2 ) for 
suitably chosen 77, fi. Here fi « |?y 2 but the most suitable 
fj, is likely to vary for other r\. The Re* m above are also con- 
sistent with the results of Laure et ah (2000) who found 
Re* n — O(10 2 ) in their calculations with fi — 0. 

4. Self-consistent dynamo 

In this section Eqs. (0) and (||) are solved simultane- 
ously using as initial conditions the previously prescribed 
Taylor-vortex flow V and the accompanying eigenfunc- 
tion B. Figure ^ shows that after an initial transient the 
magnetic energy Eb saturates to a constant value, indi- 
cating dynamo action. At the parameters in Fig. ^ the 
calculations required 16 Chebyshev modes radially, 16 ax- 
ial and 12 azimuthal Fourier modes with a timestep of 
10 -4 t, where r is the rotation period of the inner cylin- 
der. The magnetic energy is plotted as a fraction of the 
energy of the driving circular-Couette flow Eccf, the en- 
ergy source for the dynamo. The energy sink is an in- 
creased viscous dissipation in the velocity disturbance in 
addition to Ohmic dissipation. In these calculations the 
dynamo is dynamically self consistent (V and B evolve 
together). 

Figures |^ and || show the typical field structure of the 
initial conditions. A fixed outer cylinder was taken for vi- 




Fig. 8. Isosurface of helicity \V ■ V A V\ at 77 = 
Re\ = 1.5i?ei c , /i = 0. Shown over two axial periods. 



0.5, 




Fig. 9. Isosurface of \B\. The field is m 
rameters as Fig. || with P m — 2 (as in Fig. 



1. Same pa- 



sualisation purposes as the surfaces are less self-obscuring. 
The initial flow pattern V has m — symmetry and as 
Rei > Rei c vortex cores are slightly shifted towards the 
outflow regions. The eigenfunction B in Fig. ^| has m = 1 
symmetry. The flow pattern, initially axisymmetric is de- 
formed by the action of the Lorentz force (VABlAB and 
acquires an m = 2 contribution visible in Fig. |l0|. Figure 
|Tl| shows the kinetic energy of the velocity disturbance 
and the magnetic energy of the various azimuthal modes, 
E(m) and Esim) respectively, for the saturated dynamo 
of Fig. 0a. The velocity has contributions m = 0, 2, 4, . . . 
and the magnetic field has m = 1, 3, 5, . . . etc. The per- 
turbation to the magnetic field is difficult to appreciate 
visually on the dominant m = 1 structure. It remains 
rather similar to Fig. ||. 
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Fig. 10. Isosurface of helicity at magnetic field saturation. 
The flow is m = 2, looking the same if rotated by 180°. 
Parameters as in Fig. [?]a. 




2 4 6 8 10 12 

Fig. 11. Kinetic energy of the velocity disturbance and 
the magnetic energy of the various azimuthal modes. 
Parameters as in Fig. Qa. 



rated dynamo. Hopefully these results will stimulate ex- 
perimental work on the problem. Future theoretical work 
will investigate dynamo action in hydrodynamically stable 
flows and address the nature of the magnetic field struc- 
ture when the dynamo is driven harder - our dynamo is 
laminar. Most of the present work is concerned with wider 
gaps. 
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5. Discussion 

By solving the kinematic dynamo problem, we have de- 
termined that a Taylor-vortex flow pattern can sustain a 
growing magnetic field. 



Like in the models of Dudley fc James (1989) we also 



find that the dynamo is sensitive to the flow pattern. 
Further, for flows that are capable dynamo action we 
see that the growth rate is not a monotonic increasing 
function of the Reynolds number. This is not seen in 



Dudley fc James (1989) , most likely due to the prescribed 
form for the driving flow patterns. 

In the Taylor-vortex flow the best growth rates have 
been obtained with co-rotation. The relative magnitude of 
the shear and roll in the flow plays an important part in 
the success of the dynamo mechanism. 

Solving the full MHD equations we have demonstrated 
the existence of a fully self-consistent nonlincarly satu- 



